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Macroscopic price evolution models are commonly used for investment strategies. There are 
first promising achievements in defining microscopic agent based models for the same purpose. 
Microscopic models allow a deeper understanding of mechanisms in the market than the purely 
phenomenological macroscopic models, and thus bear the chance for better models for market regu- 
lation. We exemplify this strategy in a case study, deducing a macroscopic Langevin equation from 
a microscopic spin market model closely related to the Ising model. The interplay of the microscopic 
and the macroscopic view allows for a better understanding of the microscopic model, as well, and 
may guide the construction of agent based market models as basis of macroscopic price models. 
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While extensively used textbook models of financial 
markets assume an efficient market [jj and Gaussian dis- 
tributed returns of financial assets [2, S|, real markets 
show a completely different behavior with huge bubbles, 
crashes, and erratic asset price time series. This dis- 
crepancy is believed to be one of the central reasons for 
the 2008 financial market crisis Q. Scientific investiga- 
tions beyond the standard theories of financial markets 
are partly driven by methods from theoretical physics, 
especially in the description and modeling of macroeco- 
nomic quantities (e.g. price time series) as well as in the 
creation of microeconomic models with many interacting 
agents 5]. 

Price time series of markets often exhibit universal 
properties, known as stylized facts, that do not occur 
in equilibrium models of financial markets: Returns are 
power law distributed @ and their signs are uncorre- 
lated, volatility occurs clustered [7], and mnlti- fractal 
properties occur Q. The first three stylized facts have 
been incorporated in different volatility forecasting mod- 
els (broadly used for investment strategies) like ARMA 
and GARCH variants (for an overview and comparison 
see [l]). A recent development also includes multi- fractal 
properties with reasonable success [HI- Microeco- 
nomic studies have lead to a number of agent models (see 
[l^ for an overview) which, however, are mostly in the 
stage of anecdotal models. As an early model, the Lux 
Marchesi model introduces chartist agents competing 
with fundamentalists, leading to power law distributed 
returns as observed in real markets, contradicting the 
popular efficient market hypothesis. A minimalistic spin 
model 13] introduces spatial structure and exhibits pro- 
nounced bull and bear markets with a long memory as 
in real markets 
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As proposed by Farmer Q and exhibited in [T^ [13] , 
agent based models can be brought so close to real mar- 
kets that they are candidates for investment strategy 
tools. While in [l6l.[T7j the outcome of agent based mod- 
els was directly compared to price time series, we here de- 
duce a macroscopic price equation from the microscopic 
market model of ref. 141 which is based on a modified 



cedure: The focus shift from macro to micro-models is 
facilitated, the interplay of microscopic and macroscopic 
views assists agent based model construction, and finally 
with multi-fractal properties the micro- and the macro- 
representation might share the same principles. 

We here focus on the spin model defined by Bornholdt 
since it allows for a clear geometric insight, con- 
nects to the very well known Ising model, and has been 
studied as a minimalistic agent based model exhibiting 
stylized facts [l^, [l^ . We first analyze the dynamics of 
the model with a frozen coupling to the global magneti- 
zation of the system and find a phase transition between 
small and large magnetizations with a jump in volatility. 
Such a phase transition is found in other socio-economic 
models, as well "^S*]. Then we construct a macroscopic 
equation for the price, which includes a volatility switch- 
ing due to the phase transition, reminiscent of regime 
switching models The macroscopic model is com- 
pared to the microscopic one with very good agreement 
and thus provides insight in the mechanism ruling the 
spin market model. From the macroscopic point of view 
we investigate the role of the model parameters. Finally, 
we propose the use of deduced macro-equations as a guid- 
ing principle for the construction of agent based models. 

The spin model of ref. [l^ describes N = agents 
on a square lattice with periodic boundary conditions. 
The agents can decide between two states S'i = ±1 (buy 
or sell), according to two conflicting forces: On the one 
hand they tend to do what their neighbors do (as in the 
Ising model), while on the other hand they tend to join 
the minority, if minority and majority sizes diverge (as in 
the minority game [l^l). This is implemented by a local 
field 
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Ising spin model. Macro-relations for agent based models 
are discussed in [l8l - [20j and suggested for the spin model 
[3] in [2l|. There are several advantages of this pro- 



with nn(i) denoting the nearest neighbors of agent i. In 
random sequential update, a randomly chosen agent i 
evolves according to the local field hi. Its new state Si is 
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determined by the probabilities 
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For a = this corresponds to a simulation of the pure 
Ising model with a heat bath algorithm and coupling 
J — 1. If a ^ 0, the system couples to its mean field h 
as a time-dependent global quantity (comparable to the 
well known global property "price" in economy). Note 
that this coupling is of an entirely different type than 
the coupling of a spin to an external magnetic field, the 
latter being described by a term —h instead of — • Si. 
This results in a much more complex behavior compared 
to the Ising model. The time dependence of h introduces 
a feedback to the system's dynamics. Later on we will fix 
h to certain constant values in order to understand the 
role of h{t). 

The dynamics of this model exhibits large and broadly 
distributed fluctuations in the magnetization m that we 
here relate to fiuctuations in financial markets. For this 
purpose, price has been introduced via additional funda- 
mentalists, leading to an expression for logarithmic re- 
turns, which is identified with the change of magnetiza- 
tion in the model, ret(t) = \n{p{t)/p{t — 1)) oc Am [T5| . 
As can be seen in Fig. |4] on the left, where the absolute 
return distribution of the model is shown with symbols, 
the returns are distributed with a truncated power law 
and an exponent of approximately 3 (or exponent 2 for 
the cumulative distribution). In addition, a long-range 
autocorrelation as well as multi-fractal properties are re- 
ported [il]. 

In order to characterize the lattice state, the border 
line length 
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is introduced, where (ij) are pairs of nearest neighbors 
and 9(1) = 1, 6(— 1) = 0. This property is related to 
the number of agents, which are exposed to stress due 
to different neighbors and hence switch more often (The 
other agents switch with a probability close to zero. On a 
side note, this allows for a faster algorithm scaling with L 
instead of L^.). Most of the time it is near 2L (compare 
Fig. [S]) indicating that the system is undercritical most 
of the time (compare Fig. [U upper left). 

To simplify the dynamics, let us now set the mean field 
h in eq. (IT|) to a constant value (i.e. treat it as an external 
constant field) . Because the dynamics of the spin model 
leads to states with rectangular components most of the 



time, we use Si = 1 for x < ^ and Si = 



-1 for X > "I as 



the initial state. Figure [T] shows the state of the lattice 
(for the fixed time t = 400 after equilibration) and some 
time-dependent quantities for a small (left column) and 
a large h (right column). In this paper, the temperature 
is set to 0.2 Tc with Tc being the critical temperature of 
the Ising model. 
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FIG. 1. From top to bottom: Lattice states (L — 128), mag- 
netization m, and borderline length time series for h = 1.5 
(left column) and h = 2.0 (right column). 



In the case with a small constant value of h (typical 
in this model), the returns Am are Gaussian distributed 
and not correlated (also higher moments of Am are not 
correlated). Therefore we conclude that the magnetiza- 
tion performs a random walk. This can be modeled by 



Am(i) = (t(L,/i) -e* 
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with the normally distributed and independent random 
variable Eq- Q is a macroscopic equation, since it 
deals with the global quantity m, only. 

For larger h, we get a completely different behavior, 
indicating a phase transition. The system is disordered 
globally and Ih is much larger (scaling with instead 
of L). Am is much larger, as well, and it is strongly 
anti-correlated since m fluctuates fast. Fig. [2] shows the 
volatility cr as a function of h for different L, calculated 
as the average over 10000 time steps after equilibration. 
Additional boundary conditions were used in order to 
avoid states with m = ±1: All agents at a; = L/4 are 
fixed to S", = 1 (at x = 3L/4 to S^ = -1). On the left, 

the data is multiplied with (L/128)'^^^. For the ordered 
phase (small h), then the curves collapse for different L 
indicating a finite size scaling with L~^/^. On the right 
hand side, the disordered part of the curves collapses with 
a scaling cx . The phase transition is close to /icrit = 
1.86 and becomes more pronounced with increasing L. 

The reason for the scaling properties can be under- 
stood as follows: The number of active agents (those with 
much higher switching probabilities) is proportional to 
Zb. So Zb can be seen as the number of steps contributing 
to a single return, and thus the total change expressed 
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FIG. 2. Averaged volatility a for fixed values of the field h 
below and above the phase transition. L — 128 (red), L = 
1024 (blue) and L = 64, 256, 512 (grey). 
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in steps scales with , if the single steps are treated 
as independent. Including the pre- factor of in the 
magnetization and with the scaling properties of Ih (com- 
pare Fig. [1]) the observed scaling is found. 

Since in the spin market model the system is in the or- 
dered phase most of the time, a macroscopic description 
with a generalization of eq. Q is promising. The only 
change to do is to let h be time dependent. The time de- 
pendence should be smooth, so let us use a smoothened 
field h{t) = afh{t) and obtain 



Am{t) = a {L, afh{t)) ■ 
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as an equation describing the price evolution as a macro- 
scopic quantity. The dependence a{L,h) is taken from 
fits to the curves in Fig. [21 Eq. ([S]) might be compared to 
macroscopic models derived from market data [9] , which 
propose mechanisms for the generation of a{t) time se- 
ries with volatility clustering to reproduce real market's 
behavior. Here, via the frozen mean field, the macro- 
scopic equation is deduced from a microscopic model as 
ct(L, h — aifi). It is worth noting that macroscopic mod- 
els are usually constructed without using microscopic 
first principles, just with the goal of a good tunabil- 
ity and forecasting power. Multi-fractal models, on the 
other hand, include insights into the market process itself 
[lOl [T]| and thus can serve as examples for agent based 
model design, as well. 

Eq. ([S]) is valid only for fields h — am beyond the 
phase transition, because in the disordered phase anti- 
correlation takes place and the equilibration process is 
not taken into account. We can put it as follows: The sys- 
tem's magnetization performs a random walk described 
by eq. ([5]) and additionally experiences the effect of com- 
plex and finally reflecting borders around merit — ^crit/o; 
as can be seen in Fig. [3] on top. Since here we will not 
describe the action of the borders at merit, we will use 
eq. (O for supercritical m too and use the time series h{t) 
produced by the microscopic model for the further inves- 
tigation of the macroscopic model. The jump in volatil- 
ity due to the phase transition has an analogue in regime 
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FIG. 3. From top to bottom: Time series of magnetization 
m within the borders of ztrricrit; returns Am; border length 
and minimal border length allowing for m ~ 0; volatility cal- 
culated as d-{t) (black) and as a{L,arh{t)) (red). The param- 
eters are set to L = 128 and a — 20. 



switching models Q . The bottom panel of Fig. [3] shows 
the comparison of the quantities a{t) as the volatility of 
the microscopic model and a{L,h{t)) as the according 
property assigned to the macroscopic eq. ([5]) . In this pa- 
per all averages are performed over 30 time steps. The 
reasonable agreement allows for the conclusion that the 
system behaves mostly as an equilibrated system with 
the actual h and the use of eq. (O is justifled. 

Therefore, let us continue by investigating the prop- 
erties of the microscopic model using the macroscopic 
model. Firstly we discuss the distribution of absolute 
returns. Since due to eq. ([5]) a single absolute return 
I Am I > is distributed as 
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the absolute return distribution of the whole time series 
is calculated as the average 
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In Fig. |4]the absolute return density is compared to the 
outcome of eq. (jS]) with reasonable agreement. So the 
power law distribution for the case L — 128, a — 20 
emerges due to the slowly but strongly varying a{t) as a 
superposition of Gaussian distributed returns with vary- 
ing standard deviation for different times. For the larger 
system with L = 1024 and a = 60 a double humped 
structure can be seen. The additionally plotted Gaus- 
sian with the typical volatility for the disordered phase 
emphasizes its origin in the increasing difference between 
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FIG. 4. Comparison of the absolute return densities p(|Amj) 
provided by the microscopic (crosses) and the macroscopic 
model (lines). Left: L = 128, q = 20 and a power law with 
exponent —3.1 (dashed line). Right: L = 1024, a = 60 and 
a Gaussian with a — 0.0012 (dashed line). 



TABLE I. Cutoff in the density of return to zero times. 
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32 


128 


128 


512 
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10 


10 


20 


80 




328 
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Tcutofi' 


400 


30000 


6000 


20000 



returns in the ordered and the disordered phase with in- 
creasing system size. The differences at intermediate re- 
turns rely on equilibration effects around the phase tran- 
sition, which are only present in the microscopic model. 
So the mechanism producing power laws is understood 
including its limitations for larger systems. 

Finally we use the macroscopic view to discuss the 
time scales of the system. Since the partition of the ran- 
dom walk in magnetization broadens with time propor- 



tional to Am • t^l'^^ it eventually will arrive at the bor- 
ders ±/icrit/a- This defines a time scale t oc (aAra)~^ c>c 

jo? . This is in fact found in the distribution of return- 
to-zcro-times, which describes the bull and bear market 
durations. The distribution shows a power law with an 
exponent about 1/2, as can be found in the case of uncor- 
related random walks as well. The cutoff of this power 
law is shown in Table U for different system parameters 
with the finding Tcutoff oc iP" jc?, which is exactly the 
scaling we expected for the system's timescale with the 
macroscopic model. 

By characterizing the system state of the spin market 
model with the borderline length and by freezing the cou- 
pling to the global field, we found a phase transition be- 
tween small and large magnetizations with different finite 
size scaling in the ordered and disordered phases. The 
scaling was understood with the border line length. We 
deduced a macroscopic equation for the price evolution 
and compared it to the microscopic model with reason- 
able agreement, which can be understood as follows: The 
system is equilibrated almost every time according to its 
smoothened mean field. So we proceeded by investigating 
the model dynamics with the macroscopic point of view. 
With this method we found the mechanism leading to 
power laws in the absolute returns, we found the reason 
for the appearance of a double humped structure in larger 
systems and we found and explained the system's time 
scale in terms of the bull and bear market durations. In 
summary, the mechanisms driving the spin market model 
are understood in a macroscopic way and the system pa- 
rameters are characterized simultaneously. The model 
can be further adapted without losing the macroscopic 
description. So exceeding this paper, new model variants 
might be created in a controlled manner, e.g. by vary- 
ing scaling properties or including multi-fractal elements 
[13, mi that real market mechanisms suggest. 
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